Mon. Not. R. Astron. Soc. 000, [THg] (2008) Printed 14 February 2008 (MN M^X style file v2.2) 



On the relation between radial alignment of dark matter 
subhalos and host mass in cosmological simulations 



Alexander Knebe 1 , Nadya Draganova 1 , Chris Power 2 , Gustavo Yepes 3 , 

O 
O 

<N 



43 

6 



> 

On 

(N 
O 
OO 

o 



X 



Yehuda Hoffman 4 , Stefan Gottlober 1 , and Brad K. Gibson 5 

1 Astrophysical Institute Potsdam, An der Sternwarte 16, Germany 

2 Department of Physics & Astronomy, University of Leicester, University Road, Leicester LEI 7RH, UK 

3 Grupo de Astrofisica, Universidad Autonoma de Madrid, Madrid E-28049, Spain 

4 Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel 

5 Centre for Astrophysics, University of Central Lancashire, Preston PR1 2HE, UK 



Submitted Version . 



ABSTRACT 

We explore the dependence of the radial alignment of subhalos on the mass of the host 
halo they orbit in. As the effect is seen on a broad range of scales including massive 
clusters as well as galactic systems it only appears natural to explore this phenomenon 
by means of cosmological simulations covering the same range in masses. We have 25 
well resolved host dark matter halos at our disposal ranging from 10 15 /i _1 Mq down 
to 1O 12 /i _1 M0 each consisting of order of a couple of million particles within the virial 
radius. We observe that subhalos tend to be more spherical than isolated objects. Both 
the distributions of sphericity and triaxiality of subhalos are Gaussian distributed with 
peak values of (s) w 0.80 and (T) w 0.56, irrespective of host mass. Interestingly we 
note that the radial alignment is independent of host halo mass and the distribution of 
cos 6* (i.e. the angle between the major axis E a of each subhalo and the radius vector 
of the subhalo in the reference frame of the host) is well fitted by a simple power law 
P(cos9) cx cos 4 9 with the same fitting parameters for all host halos. 

Key words: galaxies: evolution - galaxies: haloes - cosmology: theory - cosmology: 
dark matter - methods: iV-body simulations 



1 INTRODUCTION 

The concordance of multitude of recent cosmological stud- 
ies has demonstrated that we appear to live in a spa- 
tially flat . A-dominated cold dark matter (ACDM) uni- 
verse (cf. ISpergel et all2007l ). During the past decade sim- 
ulation codes and computer hardware have advanced to 
such a stage where it has been possible to resolve in de- 
tail the formation of dark matter halos and their subhalo 
popu lations in a cosmological context (e.g. iKlypin et al.l 
1 19991 ). These results, coupled with the simultaneous increase 
in observational d ata (e.g . 2 deg ree Field galaxy redshift 
survey (2dF GRS). ICollessI ll200aj : Sloan Digital Sky Sur- 
vey (SDSS), lAdelman-McCarthvl <|2007f >). has opened up 
a whole new window on the concordance cosmogony in 
the field that has become known a s "near-field cosmology" 
( Free man fc Bland- HawthorrJl2003 ). 

One property of the galaxy population that has been 
measured is the spatial distribution of satellites about their 
primaries. Some evidence suggests that this is anisotropic, 
with satellites clustering about the major axis of the pri- 
mary. An anisotropic spatial distribution has been measured 



in various observational data sets (e.g. 
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2004 iBrainerdl 12005 


Yane et al.l 


2006; 


Bailin et al.l 2007 
tinues to be a m 
is obtained for tl 
of cosmological d 
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atter of debate thoug 
le distribution of sub) 
ark matter halos (e.g. 
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, as w ell as follow- 
ing from theoretical modelling ()Lee et al.ll2005l ). 

Another property of the satellite population that has 
been measured is the radial alignment of their primary 
axes with respect to the host. The first evidence for this 
effect was reported for the Coma cluster, where it was 
observed that the projected major axes of cluster mem- 
bers pr eferentially align with the direction to th e cluster 
centre l|Hawlev fc Peebles! 1 19751 ; iThompsonl Il976h . Such a 
correlation between satellite elongation and radius vector 
has further been confirmed by statistical analysis of the 
SPSS data (iPereira fc Kuhnl [SooBI; lAgustsson fc Brainerdl 
120061 ; iFaltenbacher et al.ll2007l;IWang et alj|2007ft . How erai\ 
we also acknowledge that iBernstein fc Norbera (|2002h did 
not find such a signal in the 2dFGRS. The radial align- 
ment of subhalo shapes towards the centre of their host has 
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Table 1. Details about the hosts and their subhalo populations. 
^Vp.host gives the number of particles in the host while M vlI mea- 
sures its mass (in 10 14 h~ 'Mq) and the mass -Rhost its virial radii 
(in h~ 1 kpc). The last two columns give the number of satellites in 
excess of N p > 200 and N p > 200, b/a < 0.9, respectively. These 
(sub-)samples of all identified subhalos comply with the criteria 
to reliably measure a) shape (S) and b) radial alignment (RA). 
Please refer to the text for further details. 
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also been measured for the subhalo population in cosmo- 
logical simulations l|Kuhlen et all 120071 ; iFaltenbacher et al.l 



l2007l : IPereira et alJl2007lT We no te that this effect was pre- 
dicted by Ciotti fc Duttal (1 19941 ). who used simulations to 
argue that cluster galaxies were influenced by the tidal field 
of the host cluster tidal field. 

In this Letter we provide evidence that the radial align- 
ment of subhalos in cosmological simulations does not de- 
pend on the mass of their host dark matter halo. In addition, 
we show that subhalo triaxialities and shapes follow a Gaus- 
sian distribution, which again does not depend on the mass 
of the host. 



200lf) and their properti es ha ve been discuss e d in d etail in 
Warnick fc Knebei (|2006l ) and IWarnick et all (l200St). MW1 
corresponds to model "Box20" of IPrada et al. I J2006h that 
has been simulated with the ART code (jKravtsov et al.lll997l ). 
All other halos (i.e. MC1-3, C9-12, Gl-5, M W2-5) have been 
simulated with the Tree-PM code GADGET2 l|Springelll2005h 
and their detailed properties will be presented in a forthcom- 
ing paper (Knollmann et al., in preparation). However, we 
wish to highlight that G1-G5 are drawn from a constrained 
realisation of the Local Universe (1024 3 particles in a box 
of side 64/i _1 Mpc; Yepes et al., in preparation) Q MC1-3 
are resimulations of individual halos embedded in a box of 
side 512/i _1 Mpc; MW2-4 are resimulations of individual ha- 
los embedded in a box of side 50/i _1 Mpc, while MW5 is a 
resimulation embedded in a box of side 150/i _1 Mpc. 

Note that MC1-3, C9-12 an d Gl-5 have been sim ulated 
assuming a WMAP3 cosmology l|Spergel et al.ll2007l ). while 
all other systems were simu lated assuming a WMAP1 cos- 
mology (|Spergel et alj|2003l ). 



2.2 The Subhalos 

Both host halos and subhalos are identified using 
(Knollmann et al., in preparation), an MPI parallelised mod - 
ification of the MHI0 algorithm presented in lGill et all l)2004h . 
AHF utilises an adaptive grid hierarchy to locate halos (sub- 
halos) as peaks in an adaptively smoothed density field. Lo- 
cal potential minima are computed for each peak and the set 
of particles that are gravitationally bound to the peak are 
returned. If the peak contains in excess of 20 particles, then 
it is considered a halo (subhalo) and retained for further 
analysis. 

For each halo (subhalo) we calculate a suite of canon- 
ical properties from particles within the virial (truncation) 
radius. We define the virial radius R v i T as the point at which 
the density profile (measured in terms of the cosmological 
background density pt) drops below the virial overdensity 
A vir , i.e. M{< J R vir )/(47rfl? ir /3) = A vir p 6 This prescrip- 
tion is not appropriate for subhalos in the dense environs of 
their host halo, where the local density exceeds A v j r pf,, and 
so the density profile will show a characteristic upturn at a 
radius R < R v i T . In this case we use the radius at which the 
density profile shows this upturn to define the truncation 
radius for t he subhalo . Furth er details of this approach can 
be found in lGill et all (|2004T l. 



2 THE DATA 

2.1 The Host Halos 

We use a set of 25 high-resolution cosmological (zoom) 
simulations of individual dark matter host halos. Three 
of these runs (MC1-3) represent massive galaxy cluster of 
M « lO 15 /i _1 M . We have twelve cluster-sized objects (Cl- 
12) with M ~ 1O 14 /i _1 M0, there are five group-sized sys- 
tems (Gl-5) with M w lO 13 /i _1 M , and another five Milky 
Way-type halos (MW1-5) with M ~ lO 12 ft _1 M . The par- 
ticulars of these hosts are summarised in Table [T] 

For the interested reader, Cl-8 have be en generated 
with the adaptive mesh refinement code MLAPM (|Knebe et al.l 



3 SUBHALO SHAPES 

A generic prediction of the CDM model is that dark matter 
halos are triaxial systems, that can be reasonably approx- 
imated as ellipsoids (e.g. iFrenk et al ] Il988l : IWarren et al.l 

1 See iKlypin et al . (2003) for a presentation of the method ap- 
plied to run constrained simulations. 

2 AMIGA's-Halo-Finder; AHF can be downloaded from 
http : //www. aip . de/People/aknebe/AMIGA. AMIGA is the successor 
to MLAPM. 

3 MLAPM's-Halo-Finder 

4 For a distinct (i.e. host) halo in a ACDM cosmology with the 
cosmological parameters that we have adopted, A v i r = 340 at 

2 = 0. 
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19921: Kasun & Evrardl 20051; Bailin & Steinmetz 


2005; 


Allaood et al.l 120061; Maccio et al.l 2007; Bett et al.l 


2007; 


Gottlober & YeDesI 120071^ . Following others (e.g. Gerhard 


1983; Bailin & Steinmetzl 2005 


; Alleood et al. 


2006; 


Pereira et al. 2007; Kuhlen et al. 


2007; Faltenbacher et al. 


2007h. we measure the shaDe of the halos bv the weighted 



moment of inertia: 



I ij = ^m k —^±. (1) 

h k 

The axis ratios b/a and c/a are the square roots of the 
eigenvalues (a > b > c) and the corresponding eigenvectors 
E a , Ef,, E c give the directions of the principal axes. 

We wish to measure the direction of a subhalo's princi- 
pal axes with respect to the centre of the host, and so it is 
important to measure reliably shap e and orientation. There- 
fore we follow iPereira et al.l (|2007) and require subhalos to 
contain at least N Ptin i n = 200 for reliable shape estimation 
in this Section, and to have an axis ratio b/a < 0.9 when 
investigating the radial alignments in Section [4] 

The numbers of subhalos in each host halo compliant 
with either of these two criteria are given in Table [T] One ob- 
serves that there exists a prominent population of subhalos 
with b/a > 0.9. Careful inspection reveals that the distribu- 
tion of P(b/a) (not shown here) is Gaussian with the peak 
at about 0.89. This means that there is a substantial num- 
ber of objects orbiting within the virial radius of the host 
which do not enter into ou r radial alignment an alysis below. 
However, as also noted bv I Pereira et al] (|2007f ). we confirm 
that the majority of these systems is spherical and hence a 
major axes cannot be robustly definedlj 

Because we are interested in the dependence of both 
the shape and the radial alignment on the mass of the host, 
we stack data accordingly. All subhalos for MC, C, G, and 
MW hosts are combined into one single plot and hence there 
appear four panels in subsequent plots, one for each host 
class. 
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Figure 1. Distribution of subhalo triaxialities and subhalo 
shapes. Only subalos with more than N p > 200 are considered. 



Table 2. Gaussian's parameters for the probability distribution 
of triaxiality and sphericity. 



host P(T) P(s) 





-^pcak 


<7 


•^pcak 


(7 


MC 


0.54 


0.21 


0.82 


0.07 


C 


0.57 


0.19 


0.81 


0.08 


G 


0.59 


0.18 


0.78 


0.07 


MW 


0.53 


0.20 


0.80 


0.08 



for this purpose; a spherical halo has c/a — > 1 while a highly 
oblate or prolate halo has c/a — *■ 0. The resulting distribu- 
tions for subhalos in excess of 200 particles can be viewed in 
Figure [1] (dashed histograms), too. The corresponding best- 
fit parameters of a Gaussian are again listed in Tabled 



3.1 Shape Measurements 

Triaxiality Using the eigenvalues a > b > c of the moment 
of inertia tensor w e calculate the triaxiality parameter (e.g. 
iFranx et ai]|l99ll) 



(2) 



T 



In Figure [T] we show the triaxiality probability distribution 
(solid histograms showing the fraction of subhalos in the 
respective bin and normalized to unity) for subhalos con- 
taining more than 200 particles; error bars assume Poisson 
errors. We further fitted a Gaussian to P(T) and list the two 
fitting parameters (width a and peak T pca k) in Tabled 

Sphericity Although triaxiality is useful for distinguishing 
prolate from oblate halos, we seek a more reliable measure 
for deviation from sphericity. We use the axis ratio 



(3) 



5 If we nevertheless include objects with b/a ^ 0.9 in the radial 
alignment analysis the signal is practically unaffected. 



3.2 Discussion 

There is an extensive litera ture on the shape of isolated/field 
dark matter halos (e.g. iFrenk et al.l 1 19881 ; IWarren et al.l 



19921; iKasun fc Evrardl 120051; iBailin fc Steinmet 



_ ... . — ietd 120051 ; 

Alleood et all I200& llvlaccio et al] 120071 : iBett et alJ 120071 ). 

These studies indicate that the average triaxiality and 
sphericity of isolated halos is (T) ~ 0.75 and (s) w 0.66, 
respectively. 

The first extension of such studies to s ubhalo popu- 
lations was present ed by Kuhlen et al. (|2007h for the "Via 
Lactea" simulation l|Diemand et al.ll2007t ) , an ultra-high res- 
olution simulation of a Milky Way-type dark matter halo. 
We do not have the resolution of the Via Lactea simula- 
tion but we do have a statistical sample of sufficiently well 
resolved systems to ex plore the shapes of s ubhalos. 

In agreement with lKuhlen et ahl l|2007h , we observe that 
subhalos are triaxial with (T) w 0.55 (cf. Table (2J), irrespec- 
tive of the host mass. We note also that the distribution 
of T as presented in Figure [1] (solid histograms) is well fit- 
ted by a Gaussian. The same holds for the distribution of 
sphericities P(s) pres ented in Figure U (d a shed histograms), 
too (cf. also Fig. 3 in iFaltenbacher et al.l (120071 ) ). However, 
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our sphericities are typically (s) ~ 0.80 and are marg inally 
larger than the on e s repo rted by iKuhlen et al.l (|2007h and 
IFaltenbacher et all [|2007h ((s) « 0.74). This likely reflects 
the d ifferent halo finding algo r ithms ; both IKuhlen et al.l 
( 2007) and IFaltenbacher et al.l (|2007h use a method de- 
rived from the friends-of-friends algorithm , wher eas we use 
a spherical halo finder AHF (cf. iGill et al] |2004) who per- 
formed an in-depth comparison of AHF, FOF and SKID.) 



4 SUBHALO ALIGNMENTS 

The principal aim of this study is to investigate whether 
or not there is 
of subhalos (i.e. 
respect to the 
their host halo, 
that the major 
ies tend to 



a dependence of the radial alignment 
the alignment of their major axis with 
centre of the host) on the mass of 
Recent observational evidence suggests 
axis (in projection) of satellite galax- 
point toward s the centre of their host" 
Pereira fc Kuhnl 120051 ; lAgustsson fc Brainerdl l200fj|; 



Yang et al.l 120061 ; IFaltenbacher et all 120071 ; IWang et ail 
20071 ). It is therefore natural to ask whether or not subhalos 
in cosmological simulations display a similar t rend. To date 
a few studies have investigated this subject (IKuhlen et alj 



120071 : IFaltenbacher et al.ll2007l : iPereira et al.ll2007n . We ex- 
tend this work by examining the mass dependence of the 
radial alignment. 



4.1 Measurement of Radial Alignment 

To measure the radial alignment of subhalos, we use the 
eigenvector E a which corresponds to the direction of the 
major axis a of the subhalo. We quantify the radial align- 
ment of subhalos as the angle between the major axis E a 
of each subhalo and the radius vector of the subhalo in the 
reference frame of the host: 



(4) 



R sub * E„ sub 

Rsub||E ajSu b| 

The (normalized) distribution P(cos(9) of cosO measuring 
the fraction of subhaloes in the respective bin can be viewed 
in Figure [2] 

We find a positive radial alig nment signal d ifferen t 
from isotropy, in agreem ent with IKuhlen et al.l l|2007h . 
IFaltenbacher et all (|2007l ) and IPereira et alj l|2007h . This 
can be verified in Table [3] where we compare the cumula- 
tive probability distributions P(< cos 9) (shown as thin solid 
lines in Figure [2j too) with the isotropic distribution by ap- 
plying a Kolmogorov- Smirnoff (KS) test. The resulting KS 
probabilities are consistent with zero. They are listed in Ta- 
ble |3] alongside D, the maximum distance of the actual and 
isotropic distribution. 

To better quantify these differences, we fit the following 
heuristically determined function to the (differential) prob- 
ability distributions 



P{x) = 



B + A/5 



Ax* + B 



(•>) 



with x = cos 6. In order to gauge the error introduced by 
binning the data we perform the fit to differently binned 
P(cos0) changing TVbms steadily from 5 to 15; (A) and (B) 




Figure 2. Distribution of subhalo radial alignment. The his- 
tograms are the differential distribution (with Poissonian error 
bars) that have been fitted by Eq. J5]l (dashed line). The con- 
tinuous line represents the cumulative probability distribution 
P(<cos0); the dotted line is the (cumulative) isotropic distribu- 
tion. Only subhalos with N p > 200 and b/a < 0.9 are considered. 



Table 3. Kolmogorov-Smirnov test and best fit parameters for 
radial alignment probability distribution. The errors for A and B 
are based upon fitting Eq. J5} to P(cos6>) using different number 
of bins iV bins 6 [5, 15] and (A) and (B) are the respective mean 
values. 
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KS probability 
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MC 


0.266 


6 x 10" 09 


2.81 ± 0.22 


0.56 ± 0.05 
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0.262 


6 x 10" 13 


2.67 ± 0.23 


0.60 ± 0.05 
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0.306 


5 x 10- 09 


2.50 ± 0.25 


0.60 ± 0.04 


MW 


0.327 


3 x 10" 10 


2.58 ± 0.26 


0.58 ± 0.04 



given in Table |3] are the corresponding means of the best- 
fit parameters and the errors are the lcr-deviation. We note 
that these values are practically identical across hosts and 
do not depend on mass. This suggests that there is no (hid- 
den or obvious) relation between the radial alignment of the 
subhalos and the mass of the host. This implies that any ex- 
planation of this phenomenon has to apply to galactic halos 
as well as clusters of mass 10 h~ Mq. 

4.2 Discussion 



Both I Yang et al.l (|200fJ ) and I Wang et all (|2007l ) report that 
the spatial anisotropy of satellite galaxies depends on host 
halo mass. In contrast, no such mass dependence has 
been reported for the radial alignment of satellite galax- 
ies with respect to the host. This can be understood if 
one considers the physical origin of these effects. The spa- 
tial anisotropy is linked to the anisotropic infall of subha- 
los (a nd presumably sat ellite galaxies) onto their host halo 
(e.g. iKnebe et al.ll2004h , and so can be considered an en- 
vironmental effect . In contrast, the r adial alignment is a 
dynamical effect j|Kuhlen et all 120071 ; IPereira et~ai1 120071), 
driven by the tidal field of the host halo ( Ciotti fc Duttal 
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1 19941 ); the subhalo adjusts its orientation in response to the 
tidal field on a timescale that is much shorter than the Hub- 
ble time. However, it is difficult to determine whether or 
not a subhalo adjusts its orientation by rigid body rotation 
(e.g. iKuhlen et al]|2007l ; iPereira et atl 120071 ') or by chang- 
ing its shape; it is not straightforward to distinguish be- 
tween changes in subhalo shape and "simple" figure rota- 
tion, especially for subhalos which are suffering mass loss. 
Whi le the latter has been co nfirmed for isolated/host ha- 
los dBailin fc Steinmet j [2004) it has yet to be verified for 
subhalos. 



5 SUMMARY AND CONCLUSIONS 

We have examined the whether or not the radial alignment 
of the major axes of subhalos with respect to the centre 
of their host dark matter halo is dependent on host halo 
mass. Our results draw upon a sample of 25 cosmological 
high resolution resimulations of halos spanning galaxy- to 
cluster-mass scales (lO 12 fc -1 M to 1O 15 /i _1 M ), the major- 
ity of them containing in excess of 10 6 particles within the 
virial radius. Our main results may be summarised as fol- 
lows: 

• Subhalos tend to be more spherical than isolated halos, 
with (a) « 0.80 . 

• Subhalos have triaxialities (T) ps 0.53, lower than iso- 
lated halos. 

• The probability distribution of cos 6 (i.e. the angle be- 
tween the major axis E a of a subhalo and its radius vector 
in the reference frame of the host) can be well described by 
a simple power law P(cos#) oc cos 4 9 

• These results do not depend on the mass of the host 
system. 
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